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The approximate unification of gauge couplings is tlie best indirect evidence for low-energy super- 
symmetry, although it is not perfect in its simplest realizations. Given the experimental evidence 
for small non-zero neutrino masses, it is plausible to extend the MSSM with three right-handed 
neutrino chiral multiplets, with large Majorana masses below the unification scale, so that a see-saw 
mechanism can be implemented. In this extended MSSM, the unification prediction for the strong 
gauge coupling constant at Mz can be lowered by up to ~ 5%, bringing it closer to the experimental 
value at la, therefore improving significantly the accuracy of gauge coupling unification. 
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Gauge coupling unification is an important test for 
candidates to fundamental theories beyond the Standard 
Model (SM). On the one hand, gauge unification is a 
prediction of most of such theories. In particular, GUT 
models predict the unification of the gauge couplings at 
a scale ~ Mx, corresponding to the breakdown of the 
unifying group. Also, most of the string constructions 
predict unification of all the gauge couplings, including 
gravity, even in the absence of a GUT group. On the 
other hand, due to the present high-precision knowledge 
of the gauge couplings at low energy, one can extrapolate 
them to high energy through the corresponding renor- 
malization group equations (RGEs) and verify whether 
unification takes place or not. 

Of course the evolution of the couphngs (as described 
by the RGEs) depends on the type of theory under study. 
In particular, unification can depend on whether it is su- 
persymmetric or not, whether it has the minimal matter 
content or not, whether it has (large) extra dimensions 
or not, etc. Consequently, the requirement of gauge cou- 
pling unification offers the possibility of testing these the- 
ories at very high scales. In this respect, it has been 
known for a long time that non-supersymmetric gauge 
unification (with minimal matter content) does not work. 
In contrast, supersymmetric gauge unification takes place 
with remarkable accuracy at Mx ~ 2 x 10^^ GeV. This 
fact was confirmed in 1990 with the precise LEP deter- 
minations of the gauge couplings. The unification occurs 
in such a beautiful way that it is hard to believe that it 
is a mere coincidence, and it is currently considered as 
a strong hint in favor of supersymmetry. However, there 
are some problems with the current status of this issue. 
First, the conventional prediction of (weakly coupled het- 
erotic) superstrings is that the gauge couplings (including 
gravity) are unified at the string scale Mstring ^ 5x 10^'' 
GeV, which is one order of magnitude too high compared 
with Mx- During the last years many new string (or 
string-inspired) scenarios have appeared addressing this 



problem: strongly coupled strings, large extra dimen- 
sions, etc. Another problem is that the ordinary super- 
symmetric unification, however beautiful, is not perfect 
in the simplest scenarios. A convenient way to express 
the discrepancy is to start with the SU{2) x U{1)y low- 
energy couplings ai{Mz) and a2{Mz) to obtain both the 
unification scale and gauge coupling, Mx and a{Mx)- 
Then, running from Mx back to Mz, one gets a pre- 
diction for a^{Mz), to be compared to the experimental 
value. 

Suppose then that, following for example ref. ^ one 
has a MSSM spectrum as obtained from a minimal su- 
pergravity model. Starting with the experimental input 
(in the MS renormalization scheme, before conversion to 
the DR scheme) 



a-i(Mz) = 58.98 ±0.04, 
a2'^{Mz) = 29.57 ±0.03, 
Mz = 91.197 ±0.007, 



(1) 



one integrates numerically two-loop RGEs, including also 
the supersymmetric thresholds until one finally obtains 
||l| a prediction az{Mz) — 0.13, to be compared to the 
experimental value 



az{Mz) = 0.119 ±0.004 



(2) 



which represents a 3cr discrepancy [other evaluations give 
a smaller experimental uncertainty in (^) , thus leading to 
a stronger discrepancy. We choose to be rather conserva- 
tive at this point]. Several comments are in order here. 
The effect of supersymmetric thresholds introduces a de- 
pendence on the type of MSSM spectrum considered. For 
instance, the larger the supersymmetric masses, the lower 
a-i{Mz)- Allowing the squark masses to be up to 1 TeV, 
one gets aziMz^ ~ 0.127. On the other hand, working in 
the context of gauge-mediated SUSY breaking, which im- 
plies the presence of extra (messenger) fields and thus a 
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departure from the minimal matter content, one can get 
slightly lower values, a3{Mz) ^ 0.125. This alleviates the 
disagreement, but it does not solve it. One would need a 
negative correction between— 3% and —9% in a3{Mz) [or 
a correction three times smaller in a^^Mx)] in order to 
reconcile prediction and experiment. It has been argued 
that the origin of the discrepancy could be high-energy 
threshold corrections. These may be of stringy or GUT 
origin and are model-dependent. E.g. in a GUT sce- 
nario these corrections are due to the appearance at the 
GUT scale of incomplete SU{5) multiplets, most notably 
the color triplet Higgs bosons. Demanding a threshold 
correction of the correct magnitude and sign severely re- 
stricts the models, e.g. minimal SU{5) is incompatible 
with such a requirement [Q. 

In the present paper we explore the impact of mas- 
sive neutrinos in the unification of gauge couplings. Ob- 
servations of the flux of atmospheric neutrinos by Su- 
perKamiokande provide strong evidence for neutrino 
oscillations, which in turn imply that (at least two species 
of) neutrinos must be massive. Preliminary results from 
the long-baseline experiment K2K tend to confirm this 
evidence Q]. Additional support to non-zero neutrino 
masses is given by the need of neutrino oscillations to ex- 
plain the solar neutrino flux deficit . The simplest and 
most beautiful mechanism to account for small neutrino 
masses in a natural way is probably the see-saw mecha- 
nism Its supersymmetric version has superpotential 



W = WmSSM ~ T^^RMv'li + VR^yL ■ H2, 



(3) 



where Wmssm is the superpotential of the MSSM. The 
extra terms involve, beside the usual three lepton dou- 
blets, L, three additional neutrino chiral fields vr singlets 
under the SM group (generation indices are suppressed) . 
Yy is the matrix of neutrino Yukawa couplings and H2 
is the hypercharge -1-1/2 Higgs doublet. The Dirac mass 
matrix is mo — Yi,w sin f3. Finally, is a 3 x 3 Ma- 
jorana mass matrix which does not break the SM gauge 
symmetry. It is natural to assume that the overall scale 
of A4, which we will denote by M, is much larger than 
the electroweak scale or any soft mass. Below M the 
heavy neutrino fields can be integrated out, giving rise 
to an effective mass term for the left-handed neutrinos, 
— ^v'^M^v, with Ml, = ini3^M~^mi3, suppressed with 
respect to the typical fermion masses by the inverse power 
of the large scale M. In fact, large values of the neutrino 
Yukawa couplings are perfectly consistent with tiny neu- 
trino masses for values of M sufficiently close to Mx ■ 

The influence of heavy right-handed neutrinos on the 
unification of gauge couplings is due to the fact that 
above M, the Yukawa couplings of eq.(^) affect the 
RGBs of the gauge couplings at two-loop order, in a sim- 
ilar way as the top- Yukawa coupling does in the MSSM. 
Since Y^ can be sizeable, one expects a non-trivial im- 
pact in gauge unification. The (two-loop) RGEs at scales 
Q between M and Mx, for gi, Yt and Yi, read [with gi 
normalized as in SU{5)] 



dt 

dYt 

dt 



= -2k [bi + K {bijg^ - aiatvUa)] 



(4) 



KYt (~c,gf + 3Ht + T 



K^Yt 



C'ibi + y ) g.i 



136 
I5 



glgl + 832^3 



2„2 



glgl + ( \gl + 652) Ht + (pi + i&gl ) TrHt 



dY, 

dt 



9TrH^ - 4Ht^ - 3HtT - 9TrH. 
kY, {-c'^gl + 3H, + T) 



(5) 



c'M + ^ gt + plgi 



&\\gl + gi) H. + (\gl + 16.g| ) TrHt 



STrH^ - ml - 3H^T - 9TrH^ 



(6) 



(summation over repeated indices is understood) with 
t = InQ, K = l/(167r2) and 



H„=YtY„, T = Tr(3Ht + H,), 



(7) 



with a ^ Here Yt represents the Yukawa matrix of 
the It-type quarks (dominated by the top Yukawa cou- 
pling) and the numerical coefficients are given by 

h = (33/5, 1, -3), c, = (13/15, 3, 16/3), 

26/5 6/5 
= (3/5,3,0), a,„ = I 6 2 
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199/25 27/5 88/5 \ 
6y = I 9/5 25 24 

11/5 9 14 y 



(8) 



It is clear from eq.(|j) that the presence of neutrinos will 
only affect the RGEs of gi and g2. The magnitude of the 
final correction will depend on the numerical values of Y^ 
and M . These are not independent quantities, since the 
final neutrino masses My = v'^ s\t? (3Y M~^Y ^ (ap- 
propriately ran down to the electroweak scale ^) must 
be consistent with observations. Atmospheric neutrino 
data only allow to determine one neutrino mass splitting, 
AmJ '-^ eV^. On the other hand, 'standard' expla- 

nations of solar neutrino fiux deficits require neutrino os- 
cillations with a much smaller mass splitting, correspond- 
ing to a different pair of neutrinos. In addition, there are 
upper bounds on neutrino masses, coming e.g. from the 
non-observation of neutrinoless double /3-decay and other 
experiments. All this information implies that there are 
three possible types of neutrino spectrum: Hierarchical, 
m\ < ml ^ TO§; Intermediate, m\'^ and De- 



generate, 



In the hierarchical case, the 



larger neutrino mass should be 0(10 ^ — 10 ^) eV, while 
in the degenerate case cosmological observations require 
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nil, ~ C'(2) eV. This means, that once a particular sce- 
nario for the neutrino spectrum is chosen, Yj, and M are 
not independent quantities any more. 

The presence of right-handed neutrinos above M in- 
duces a small change both on the scale of unification and 
the value of the unified gauge coupling. This in turn 
affects the prediction of a3{AIz) derived from the as- 
sumption of perfect unification. Given that the changes 
with respect to the MSSM case can be treated as a small 
perturbation, it is simple to make an analytical estimate 
from eqs.(^ and (||) to get: 



1 <TrYtY.>ln^ 



M 



(9) 



In this result, < TrYj^Y,^ > is an average value in the 
interval from M to Mx- 



< TrYt Y^ > In 



M 



Mx 



TrYtY, d\nQ. (10) 



Q=M 



Two important implications of Eq. (|^) arc, first, that 
neutrino corrections always make the predicted a^lMz) 
smaller, as required to improve the agreement with the 
experimental value, and second, that the effect will be 
important if the neutrino Yukawa couplings Yj^ are siz- 
able, which is a natural assumption, as we saw. 

We can now estimate the average in Eq. ( |lo|) consid- 
ering that the evolution oiY^{t) is well described by its 
1-loop RGE neglecting all couplings different from Y^. 
This approximation is justified if Yi^ is large (the case of 
interest). The final result will depend on how many neu- 
trino Yukawa couplings (or more precisely, eigenvalues 
of Yi^^Yi,) are large. With only one large neutrino cou- 
pling (this corresponds naturally to a hierarchical neu- 
trino spectrum) it is straightforward to evaluate the in- 
tegral in (Wm in the approximation just described to get 



37V 

ZOTT 



Y^My 



(11) 



with A'^ = 1 counting the number of large neutrino 
Yukawa couphngs. Of course the effect on a3{Mz) is 
larger if all three neutrino Yukawa couplings are large 
(this corresponds naturally to a degenerate neutrino 
spectrum, but not necessarily: a hierarchy in the spec- 
trum can be induced also in this case by flavor-dependent 
Majorana masses). In this = 3 case, the result for 
Aa3{Mz) slightly depends on the kind of texture for Y^. 
A good estimate isgiven by simply taking Yi, =Y^1 and 
in that case, Eq. (|ll| ) also holds, simply setting = 3. 

The magnitude of the correction (^l]) depends on the 
values of Y^, at Mx and M, but the latter is fixed in or- 
der to fit the physical neutrino mass, as explained above. 
Hence, the correction is just a function of a unique pa- 
rameter, y^(Mx). [It is also a function of the neutrino 
mass, m^, but changes on just imply a corresponding 
modification of M and since the evolution of Yi, occurs 



mainly near Mx the final correction on a3{Mz), given by 
eq.(|l^), is quite insensitive to m^.] The most favorable 
case occurs when Yi,{Mx) is large, near a Landau pole. 
A representative case is Y^{Mx) ~ 8. Then Y^{M) - 2 
and we get /S.a-i{Mz) ~ —1.4% for the case = 1 and 
Aq;3(Mz) ~ -4.1% for A^ = 3. It is clear that in the first 
case this correction is not sufficient to reconcile aj,{Mz) 
with the experimental value, but it can be enough in the 
second case. These estimates agree well with our full 
numerical results, obtained by numerical integration of 
the RGEs (Q-|^). Now the two-loop contributions soften 
the RGE for Y,y, which implies that Y^, decreases more 
slowly when the energy scale goes down. This increases 
the Yukawa coupling average in eq.(^), and thus the (neg- 
ative) shift of the low-energy strong coupling, A.az{Mz)- 
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FIG. 1. Unification of the gauge coupUngs o"^ for a typical 
MSSM scenario, without massive neutrinos (upper plot) and with 
them (lower plot). 

An important question at this point is how large can 
Y^{Mx) be without jeopardizing the perturbative expan- 
sion, and hence the reliability of the results. This has 
been studied in great detail in rcf. jl^ for the case of the 
top Yukawa coupling. Since, in the large coupling regime 
relevant for this question, the RGEs for the top and neu- 
trino Yukawa coup ling s (in the A'^ = 3 case) are identical, 
the results of ref. |l^ are applicable also to the present 
case. The conclusions of the authors of ref. ||l^, us- 
ing Pade-Borel resummations of the 4-loop beta-function, 
were that the qualitative behaviour of the running is well 
described by the one-loop approximation and can be fur- 
ther improved by Pade-Borel approximants which are re- 
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liable for values of the Yukawa coupling up to ;$ 8 (For 
the = 1 case this is a conservative estimate, since in 
that case the beta-function for is softer than that of 
Yf.). We have used a [1,1] Padc-Borel resummation of 
the beta-function to obtain our most reliable numeri- 
cal results. As the beta-function in this approximation is 
smaller than at one-loop, the final result for the average 
in (^) turns out to be a bit larger than our analytical es- 
timates [for a fixed value ofY^{Mx)]- The negative shift 
on a3{Mz) goes then from 4% to about 5%. This fact, to- 
gether with the uncertainties in ai{Mz), a2{Mz) and the 
top mass, imply that the usual MSSM scenarios, which 
yielded an unsuccessful prediction a'i{Mz) — 0.13, are 
now rescued within la confidence level. [The agreement 
can be perfect in other MSSM scenarios, which predict 
slightly lower values of a'i{Mz), even if the error in (|^) 
becomes smaller.] 

This is illustrated in fig. 1 for a typical MSSM sce- 
nario The first plot shows how the gauge couplings fail 
to unify in the ordinary MSSM in the absence of mas- 
sive neutrinos (the requirement of unification would im- 
ply a^{Mz) = 0.13 in the case plotted). The second plot 
corresponds to the MSSM extended with neutrinos get- 
ting mass via a see-saw mechanism. More precisely, the 
plot corresponds to a scenario with degenerate neutrinos 
of mass = 2 eV and Yi,{Mx) = 8. It is apparent how 
the tti and a2 runnings are modified in a suitable way to 
get gauge unification. 

Besides the logarithmic effect on unification we have 
described, the presence of heavy right-handed neutrinos 
affects the running of the gauge couplings also through 
finite two-loop threshold effects at the Majorana scale 
M . However, these will depend on the neutrino Yukawa 
coupling y^(M), which in the case of interest is much 
smaller than Yi,{Mx), so that it is safe to neglect these 
matching effects. 

Finally, it is interesting to point out that effects similar 
to the ones we have found are expected in the next-to- 
minimal supersymmetric standard model (NMSSM), i.e. 
the MSSM extended with a singlet chiral multiplet S. 
The superpotential of the NMSSM does not contain a 
mass term for the Higgs multiplets (the /x-term). How- 
ever, an effective /i term of the correct order of magnitude 
is generated dynamically by AVF = YgSHi ■ H2, when S 
takes a vacuum expectation value. This solves in fact the 
/i problem of the MSSM and is one of the main virtues 
of the NMSSM (although this model has its own draw- 
backs). The influence of the new Yukawa coupling Kj on 
the running of the gauge couplings is also a two-loop ef- 
fect, of exactly the same form as in (||) with Y^ instead 
of TrHjy (the coefficients Uai are exactly the same for 
a = v and a = s). The final impact on a^{Mz) is there- 
fore given by a formula like (|) with TtYIY^-^Y^ and 
M replaced by the mass Ms of the singlet 5", which is 
close to the electroweak scale (this makes the logarithm 
much larger). Numerically, we find /^a^{Mz) — —3% for 
Ms = 1 TeV and Ys{Mx) = 8. Similar effects can occur 
in other extensions of the MSSM with additional Yukawa 



couplings, as has been shown for models with i?-parity 
violating couplings in ref. [^. 

In conclusion, we have examined the impact of heavy 
see-saw neutrinos (plausible in view of the growing exper- 
imental evidence in favor of non-zero neutrino masses) 
on the unification of gauge couplings, in particular as re- 
flected in the uniflcation prediction for the strong gauge 
coupling constant at the electroweak scale. We find that 
the effect is small, but is always of the right sign and can 
be of the right magnitude to bring the too high MSSM 
prediction for a^^Mz) down to values within la of the 
experimental value. Given that adding three heavy right- 
handed neutrinos is not an ad-hoc extension of the MSSM 
but on the contrary is well motivated by experiment and 
theory alike, this result is welcome and noticeable. This 
effect should be taken into account, even in models with 
sizeable stringy or GUT high-energy threshold correc- 
tions. For example, models that have been discarded 
for not giving the appropriate threshold corrections (e.g. 
minimal SU{5), [jl|), can be now perfectly consistent. 
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